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Abstract 

This paper is understood as a supplement to the paper by [Stutzki et 
al, 1998] , where we have shown the usefulness of the Allan- variance and its 
higher dimensional generalization, the A-variance, for the characterization 
of molecular cloud structures. In this study we present the connection 
between the Allan- and A-variance and a more popular structure analysis 
tool: the wavelet transform. We show that the Allan- and A-variances 
are the variances of wavelet transform coefficients. 



1 Introduction 



In | Stutzki et al, 199§}] (hereafter paper I) we defined the A-variance as a gen- 
eralization of the Allan- variance, a method traditionally used in the stability 
analysis of electronic devices. The Allan- or A-variance method allows a clear 
separation between the global drift characteristics of the signal and other con- 
tributions, like white noise. It is specially suited to estimate the spectral in- 
dex of power law spectral distributions. We demonstrated the use of the two- 
dimensional A-variance for the analysis of molecular cloud images, e.g. maps 
of CO line emission. Application of the A-variance analysis to the integrated 
CO intensity images of various molecular clouds shows a power law behaviour 
of power spectrum of the intensity distribution with spectral index (3 close to 
2.8. This intensity distribution has a structure well described by a fractional 
Brownian motion (fBm) model and in paper I it was shown that the power law 
index (3 derived by the A-variance is related to the drift exponent H of the 
corresponding fBm model by 8 — E + 2H (where E is the Euclidean dimension 
of the model) and to other parameters characterizing the fractal, such as the 
area-perimeter index. 

Another tool well established over the last decade in the areas of struc- 
ture analysis and data processing is the wavelet transform. While the dis- 
crete wavelet transform is mainly used in the field of signal compression, the 
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continuous wavelet transform has applications in texture analysis and feature 
extraction. The connection between the local regularity of a function and 
the scaling behaviour of its wavelet transform coefficients, as the scale con- 
sidered becomes smaller and smaller, makes the wavelet technique a natural 
tool for the analysis of self-similar distributions with (multi-)fractal properties 
(eg. [Holschneider, 1995 ). 

In the context of molecular cloud struc ture analysis the wavelet transform 
was first used by |Gill fc Hcnrikscn, 1990 , who also pointed out the advan- 
tages of this method over traditional structure indicators like the autocorrela- 
tion function. In their study they found a power law behaviour for the centroid 
velocities of 13 CO s pectral data of L1551, as well as for the peak intensities. 
[ Langer et al., 1995 ] , presented a clump decomposition of CO maps of the cloud 
Barnard 5 using the Laplacian pyramid transform, which is also a wavelet trans- 
form. The number distribution of clump masses derived by the Laplacian pyra- 
mid transform follows a power law, in agreement with the results presented by 



[Kramer et al., 1998 , obtained by using gaussian clump decomposition. These 
results give additional evidence for the hierarchical structure of molecular clouds. 

The purpose of this paper is to show that both these concepts are closely 
related, namely that the Allan- and A- variance can be expressed using the 
variance of wavelet transform coefficients for suitably chosen wavelets. 

After a short introduction of the Allan- variance and the wavelet transform in 
one dimension we show the relation between these concepts (section 2). The use 
of the A-variance and the continuous wavelet transform generalizes our result 
to two dimensions (section 3). We summarize our results in section 4. 



2 Allan-variance and wavelets 
2.1 Allan- variance 

We first repeat some facts about the Allan- variance; for a broader introduction 
the reader may consult appendix A of paper I. 

To calculate the Allan- variance of a random signal s(t) we consider the dif- 
ferences d(t,T) of succeeding averages over a (time) interval T and calculate 
their variance. Defining the down-up rectangle function by 




- T < t < 
< t < T 
elsewhere, 



this difference can be written as a convolution 




-{s*2^ T }(t) : = d(t,T). 
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Without loss of generality we here assume the average s to be zero. Then the 
Allan- variance cr\{T) is defined as 

a\{T):=\{£^T)) u (2) 

where (...) t denotes the average over all times (or positions) t. 

2.2 Wavelets 

The Fourier transform of a function ip is defined by ip(iL>) — J JR e^ l2 ' KUJt ip(t)dt. 
A square integrable function ip, which satisfies 

< := / S^cfc, < oo (3) 
J jr. M 

is called a wavelet. For a 6 M\{0},6 6 1R the wavelet transform of a signal s 
with respect to ip is defined by 

L^s(a,b) := ^lal" 1 / 2 / s(t)yj (—} dt. (4) 
V^p J JR. \ a J 



Each wavelet transform coefficient L^s(a, b) is thus the scalar product of the 
signal with a dilated and translated version of the "mother function" ip. 

2.3 The Allan- variance as variance of wavelet transform 
coefficients 

As a special case of a mother function consider a translated and dilated version 
of the well known Haar wavelet tpjj, which we define here as follows: 



<Y,(/): { -1 




-1 < f < 
< t < 1 
elsewhere 



For this the constant defined by (||) is c$ B = 8 In 2 (see Appendix). In the 
case of this modified Haar wavelet the transform coefficient L^s(a 7 b) given by 
(H) corresponds to the difference of succeeding integrals over s(t), which is the 
same as the averaging process leading to the definition of the Allan- variance. It 
follows that 

b b+a 

ipfi (j—) s< y l ) dt = J S ^) dt ~ J s(t)dt, 

JR. b-a b 

which equals the term in brackets in eq. (Q) with b = t and a = T. Further we 
have 

L^ s s(a,b) = -^-\a\- 1 ' 2 [ s (t)ip(—)dt 

J<^~ JjR V a J 
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i.e. d(6,a) = J —Z-L^, R s(a,b). (5) 

Using eq. (||) it follows that the Allan-variance can be expressed as 

al(a) = ^{Ll_.s(a,b)) b , (6) 

that is the Allan-variance at lag a equals the average of the wavelet transform 
coefficients over the translation index b for constant dilation index a. 



3 A- Variance and A- Wavelet 

The main idea of the proceeding section was, that the difference of adjacent 
averages corresponds with the convolution product of the signal s with an oscil- 
lating function tp of zero mean, that is a wavelet transform coefficient L^s. This 
concept will now be extended to the symmetric version of the Allan- variance, 
the A-variance as introduced in paper I. We recall from paper I (Appendix A, p. 
713 ff.) the definition of the one-dimensional down-up- down rectangle function 

f 1 : M<l/2 
iJUaO := \ -1/2 : 1/2 < \x\ < 3/2 . 
[ : elsewhere 

The double difference A(t,T), that is the average over succeeding intervals is 
defined as follows: 

A(*,T): = {*(...) *^JL(f)}(«) 

= / <*4JLi (^) s(x). (7) 

TR. 

The A-variance is the variance of the double difference A(t,T): 

a!(T):=i(A 2 (*,T)} t . (8) 

To find the relation to a wavelet transform we choose the down-up- down rect- 
angle function |_f~|_| (a:) as mother function, which defines the family of one- 
dimensional "A- wavelets" Ft 

*pA a ,„( x ) '■= JL (~~) ' 

1 A scaled version of this A- wavelet is used as "French Hat" wavelet in the literature, cf. 
[Gill & Hcnrikscn, 1990 1. 
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where again a 6 H/{0} and b E IR. The value of the admissibility constant is 
c^ A w 3.37. The corresponding wavelet transform coefficients are 

L^s(a,b) = ^|a|- 1/2 / dxdh (—)•*(*), 

a/C^I J V « / 

IR 

and with a = T,b = t and lTLK^) = lTLKM) one gets 



A(6,o) ( =' J^-L^s(a,b). 
With cq. (||) this leads to the result analogous to eq. (^|): 

<ri(a) = ^(Ll A s(a,b)) b . (9) 

3.1 Two-dimensional Continuous Wavelet Transform 

To generalize our results to higher dimensions we make use of the definition of 
two-dimensional (directional) wavelets given by [Antoine & Murenzi, 199E]. As 
in the one-dimensional case a two-dimensional function tp is called a wavelet, if 
the following admissibility condition holds: 

0<c»:=/ W*r<~. 
J1R 2 Il/H 2 

From [ Antoine fc Murenzi, 1996| we adopt the notation for translation, dilation 
and rotation operators: 



translation: (T b s)(x) 
dilation: (D a s)(f) 
rotation: (R e s)(x) 



s(x-b), 6elR 2 

±s(§), a>0 
s(r_ e (x)), 6 G [0, 27r), 



where r £ SO(2). A combination fl(a, 0, b) — T h D a R 9 of these operators acts 
on a function s in the following way: 

(fl(a,e,b)s)(x) = s a f) i {x) = ~s (~r- e (x - 6) J . 

With this notation the two-dimensional continuous wavelet transform (CWT) 
of a function (or "image" ) s with respect to ip is given by 

L^s(a, 9, b) := ^ / d 2 f V ( r_*(— ) ) (10) 

with a € IR\{0}, 6 S IR 2 . For a complex valued mother function tp in eq. ( |l0| ) is 
replaced by its complex conjugate ip* , but we will only make use of real valued 
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functions here. At present only isotropic wavelets ip{f) = ip{r),r = \r\ will be 
considered, so that we omit the index 8: 



L^s{a,b) := 



d 2 x ip I — — - ] s(x). 



3.2 A- Variance in Two Dimensions 

This section is based on paper I, Appendix B (p. 717), where we now consider 
the case E = 2, that is two-dimensional signals or images. 

The two-dimensional down-up- down rectangle function is given by: 



l_TL 0" 



(3 2 - 1) 



1 

-1/8 




1 r 
n (r) — ; n(-l 

\f\ < 1/2 
1/2 < \ f\ < 3/2 
elsewhere 



For the two-dimensional unit sphere the "volume" is V2 = it and the "surface" 
is 1S2 = 2ir. Using the nomenclature of paper I we have 



D,2 



(11) 



Convolving s(r) with Q D2 ( r l leads to the difference of averages at average 
distance D: 



A s (r,A2) := U...) * (f) 
= / d 2 iQ(r - x) ■ s(x) 



TR' 2 



D.2 



-(f)' 



d'xdht 1 - 



D 



s(x) 



(12) 



IR 2 



As in the one-dimensional case the A-variance is given as the variance, i.e the 
autocorrelation function at zero lag of the filtered signal: 



(13) 



3.3 A- Variance and Two-dimensional Continuous Wavelet 
Transform 

To relate the A-variance to the two-dimensional CWT one naturally considers 
the down-up- down-rectangle function given by eq. (O) as mother wavelet and 
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defines 



ipA„ ,(r) :=|_TLl 



b 



V J A - will also be called (two-dimensional) A-wavelet; the admissibility constant 
for the two-dimensional function is c^ A ss 4.11. Like in the one-dimensional 
case the wavelet transform coefficients L^, A s(a,6) are directly related to the 
differences of average values at average distance D A(r, D, 2) (cf. eq. (fiT 



L^ A s(a,b) = 



1 1 



d 2 x -0a -{x) ■ s(x) 



m 2 



1 1 



/c,/> A a 



- d 



We identify D = a and r = b and because of |_fl_j ( 
eq. © and eq. @: 




L. 0A s(a,&) 



1 1 7ra 2 . ,t> 

— — -A. M,2) 

/c^, A a 4 

A»(6,o,2). 
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Looking at the autocorrelation function as the two-dimensional generalization 
of the variance as we have done in eq. ( |i"3[ ) , we find as the main result of this 
section: 



2 i \ 1 16 c^, A 

Z7T 7T z a z 



s(a,b) 



o ,6)(0), 



(15) 



i.e. the A-variance for a given lag a can be expressed in terms of the autocor- 
relation function at zero lag of the wavelet transform coefficients with dilation 
parameter a. 



4 Summary and Outlook 



The main results of our study are equations (||), (||) and ([L5|) which show, 
that the Allan- and A-variance can be expressed using the coefficients of spe- 
cial wavelet transform s and that both these concepts are very closely related. 
[ Pando fc Fang, 1998 have shown that the discrete wavelet transform is a good 
power spectrum estimator especially in the case of finite sized samples. The 
same result was found by [Bensch et al.,m. prep.,\ for the A-variance analysis 



of CO integrated intensity maps as small as 32 x 32 pixel. Detailed investiga- 
tions by Bensch et al. and [Ossenkopf et al., in prep.\ show that the A-variance 



remains a robust measure of astronomical images if the effects of noise, smearing 
due to the finite size antenna beam pattern and radiative transfer are taken into 
account. 
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Pando ct al. pointed out that evaluating the power spectrum with an over- 
complete basis function system of the C WT can lead to correlations " that are 
not in the sample, but due to correlations among the wavelet coefficients" . This 
will also be true for the Allan- and the A-variance, because the definitions of 
both employ such an over-complete system. It has to be checked in a further 
study how the results change, if only coefficients on a discrete, wavelet adapted 
grid (a "wavelet frame") are considered, i.e. if the continuous wavelet transform 
is replaced by a discrete one. 

While the choice of isotropic wavelets was natural and sufficient in our work 
on the isotropic scaling behaviour of molecular cloud structure, the directional 
sensitivity of non-axial symmetric wavelets will allow the determination of pre- 
ferred orientations, eg. of filamentary substructures, in astronomical images. 

A Evaluation of the factors 

For all wavelets discussed in this paper the factors 



for the one- and two-dimensional case, respectively, can be calculated analyti- 
cally. The integrals were solved using Mathematica 3. for students. 




and 




A.l Modified Haar- Wavelet 



Looking at the mother function ip(x) = — 2|_P(a;) of the modified Haar wavelet 
system we find (paper I, p. 714) 



I 7/ \ 1 2 a S1U V 711 -^ I i i r 

=4—^ ttj— and therefore 





A.2 A-Wavelet 



A. 2.1 A- Wavelet in One Dimension 
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^ sin 6 (7rcj) 
(ttoj) 2 ■ \u\ 



47T 



sin 6 (7rcj) 



91n3 

12 In 2 — w 3.37 



A. 2. 2 A- Wavelet in Two Dimensions 

The power spectrum of the two-dimensional down-up-down rectangle function 
is (paper I, p. 717) 



J^irf) Ji(3tt/) 



(*/) (37T/) 

where Ji denotes the first order SesseZ-function of first kind. With this we get 



|JL(/)II 



■d 2 f = 2tt 



DC- 

/ 



Jl(*7) _ ■/! (37T/) 

W) (3tt/) 



-4f 



]R 2 

2tt 



(ln(3)-i) 



4.11 
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